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' Abstract. The quantum discrete Liouville model in the strongly cou- 

f"^ , pled regime, 1 < c < 25, is formulated as a well defined quantum me- 

£Nj ■ chanical problem with unitary evolution operator. The theory is self- 

dual: there are two exponential fields related by Hermitean conjugation, 
satisfying two discrete quantum Liouville equations, and living in mu- 
tually commuting subalgebras of the quantum algebra of observables. 
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Introduction 



The Liouville equation [33] 



(1) 4>tt - <t> xx - 4e" 2 * = 



has plenty of important applications in Mathematics and Physics. In par- 
ticular, it describes the surfaces of constant negative curvature and plays 



the indispensable role in uniformization theory of Riemannian surfaces [35 
(for some recent approach see EH, 41, 42|). 



O ■ In modern physics the Liouville equation defines one parameter family 

<D of models in conformal field theory (CFT), which usually is identified with 

2-dimensional gravity [|30[|. It plays even more important role in Polyakov's 
theory of noncritical Bosonic string in dimensions d < 26 p6]|. For these 
reasons the Liouville model, especially in its quantum version, attracted 



5_i ■ wide attention during the last 25 years |10|, |11 



The parameter c, labelling the quantum Liouville theory as a CFT model, 
is the central charge in a representation of the Virasoro algebra. The quasi- 
classical (or weak coupling) region, corresponding to large positive c, is well 
understood. The domain c < 1, containing the minimal models of CFT for 



some preferred discrete values of c, is also well described ||, 22, 18]. It is 
the region 1 < c < 25, corresponding to strong coupling, about which there 
exists very limited knowledge until now. 

In this paper we begin to describe one more method to treat the Liouville 
model, which is applicable for studying the strong coupling region. The 
method is based on the apparatus of quantum integrable models (see e. g. 



1 14 ] for the recent survey) . Some parts of this machinery were already used 
in CFT and, in particular, for the Liouville model. However, we feel that 
we have added several new things to this development. 

First, we use the lattice regularization for the model, which exactly retains 
the integrability. Here we follow the previous papers |l6]. |2l|| . Second, we 
show that it is indispensable ( especially in the strong coupling region) 
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to use simultaneously two mutually dual discrete models, corresponding to 
two exponents of the coupling constant r , symmetrically entering the 
expression for the central charge 

c = i + 6(t + t~ 1 + 2). 

Positive r correspond to c > 25 or weak couplings, while strong couplings 
lead to complex r on unit circle, 

r = e w , T - 1 = e- i6 = f, 

and only unification of the models for r and 1/r can restore unitarity. Sim- 
ilar considerations on constructing "modular doubles" were used earlier in 
simpler examples of the Weil-Heisenberg algebra and quantum group 
||. This paper, technically being more involved, ideologically belongs to 
the same line of thought. We must stress, that one can see the elements of 



the dualization idea in papers M, uA 10, devoted to the Liouville model 



itself, and in [||, ||, ||] within the framework of conformal field theory. 

Our main result consists in constructing the unitary evolution operator 
for the chiral shift serving the both dual models. 

In the first part of the paper we remind various facts about the Liouville 
equation to be used in what follows. Many of these facts are known in one 
or another form; however we present them in the form most suitable for our 
goal. Then follows description of the discretized Liouville equation and its 
formal quantization in purely algebraic manner. The appropriate involution 
together with proper quantization are introduced in section The self-dual 
structure is indispensable for that. 

This paper is partly supported by RFBR grant 99-01-00101, grant INTAS 
99-01705, and Finnish Academy. A.V. wishes to acknowledge the financial 
support extended by the DWTC office of the Belgian government through 
the IUAP project P4/08. A.V. and L.D.F. are also grateful to ESI, Vienna 
for hospitality during October, 1999. 

1. Recollecting the facts 

In this section we remind in appropriate form some facts about the Liou- 
ville model and the technique of its discretization. 

1.1. Liouville formula and Mobius-invariance. Liouville formula makes 
solutions, all of them in fact, from pairs of arbitrary functions of one variable 
"moving" in light-cone directions 

(2) _ P -2H*,t) = a'(x-t)(3'(x + t) 

1 ' {a(x - t) - 0{x + t)) 2 ' 

The right hand side of eqn (||) is invariant under simultaneous point- wise 
Mobius transformations 

m 2 i+am 11 no** 
m 2 2 + am\2 

of the "chiral halves" a and (3. Accordingly, if those are "Mobius-periodic" 
a(x + 2tt) = a(x) < T (3(x + 2tt) = f3(x) < T, 
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with the same monodromy matrix T in both chiralities of course, the solution 
comes out periodic in the spatial direction 

(j>{x + 2ir,t) = 4>(x,t). 

From now on, this will be the only boundary condition in use. 

Because of this hidden Mobius-symmetry, the chiral halves cannot be 
uniquely restored from a given solution, but their Schwarz derivatives can, 
for they are Mobius-invariant as well. In particular, that of a (times minus 
one half) 

doubles as the (chiral component of the) stress-energy tensor 

(3) u(x -t) = \{<t> x - 4> t ? + \^>xz ~ <kx) + e" 2 ^. 

We will now pay less attention the parallel, or rather perpendicular, j3- 
chirality. 

1.2. Magri bracket. The canonical Poisson bracket 

{w(x),ip(y)} = 5(x-y) 
on the laboratory Cauchy data 

<P = 4>\t=o ™ = <j>t\t=o, 

together with canonical total momentum and energy 

P = J w<p'dx H = \ J(w 2 + ip' 2 + ie~ 2lp )dx 

gives the equations of motion in the Hamiltonian formalism: 
if' = {P, ip} ip = w = {H, (p} 

m' = {P, w] w = ip" + Ae~ 2v = {H, w}. 

These equations lead to a single free motion equation 

(4) u = -v! = {L ,u}, 

with Lq denoting the (slightly shifted) zero'th Fourier coefficient 

L Q = j\u + \)dx = \{H -P + tt). 

The bracket involved is the same canonical bracket of course, but now we 
are prompted to express it, using formula (|3|) at t = 0, in terms of either u 
itself or its Fourier coefficients L a = f{u + j)e~ %ax dx. The Magri bracket 
II 

(5) {u(x),u(y)} = (u(x) + u(y)) S'(x - y) - \ S'"(x - y) 

and yet more celebrated (Poisson bracket realization of) the Virasoro algebra 

-i{L a , L b } = (a- b)L a+b + 7r(a 3 - a)S a _ 6 

emerge. They will attract our attention in the next couple of (sub)sections 
where we shall try to treat the Magri bracket and the "chiral Hamilton- 
ian" Lq as a stand-alone dynamical system, not for the first time of course. 
We shall review some earlier attempts to get it properly discretized and 
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quantized. It will give an idea how far we can go on bare common sense and 
formal algebra, without actually touching the Liouville equation. 



1.3. Volterra model. The first such attempt [16] produced a viable lattice 
counterpart of the Magri bracket 

{"m]"n} — 4 U m U n {^{4: U m Un) (<^m+l,n &m— l,n) 

- Ul( m+n) {5 m+ 2,n ~ 5 m -2,n), 

with u n ~ 1 — A 2 u(nA). In the traditional difference-differential scheme of 
things it has to do with Volterra's venerable preys and predators 

(It) u n = u n (u n -i -u n+ i) = log u m ,-u n j , 

otherwise known as the lattice KdV equation. It should be noted however 
that the latter only emerges in a rather tricky continuous limit, while the 
simplest limit leads to the free equation (|j). Which is a little unfortunate, 
for however essential a part has the Volterra system been playing in the 
soliton theory, there must be something wrong about a nonlinear equation 
emulating a linear one. Apparently, discrete space and continuous time do 
not get along, still, the KdV connection may teach us a thing or two. 

1.4. Bi-hamiltonity. Historically, the Magri bracket brought about the 
notion of "raising and lowering" which soon became paramount for the whole 
Hamiltonian theory of soliton equation. In our case, one lifts Lq to 

Lq = J \ u 2 dx 

to turn the free motion (|j) into the KdV equation 

„; - I,,'" _ in,*, 1 - /r.T 
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2>uu = {Lq, u}, 



then drops the Magri bracket to that of Zakharov-Faddeev's 

{u(x),u(y)}i = 25' (x - y) 
to restore the free motion in a different guise 



u = -v! = {LJj"};- 



Although not quite relevant to the Liouville equation, this makeshift may 
be useful anyway, for the corresponding lower guise of the Volterra system 



Un = U n (u n -\ - Un+i) = |y^U> 



features a more quantization-friendly bracket 

{^mjU-nij. — 1± m 11 n (<^n+l,n $m— l,n)- 

We happen to already know [jE]] how to turn this one into a noncommutative 
algebra and to arrange there a "free motion in discrete time". So, in the 
rest of this section we shall be looking at that motion, then we shall try to 
guess how it could be lifted back to the Magri- Virasoro case. 
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1.5. Quantization. Going quantum, bracket j_ ought to become Weyl-style 
exchange relations 

u "m u "n — y ""n^m • 

We do not want to discuss before time the nature of quantization constant q 
or the exact degree of formality implied in what follows, but we want to be 
very clear about boundary conditions: the u's are strictly periodic and so is 
the Kronecker symbol of course 



U n +N = Ur, 



1 if m = n (mod N) 
otherwise, 



on top of that, for the reason to surface shortly, we want period N to be 
even and an additional condition 

(6) U\U 3 . . . UN-l = U2U4 . . . UN 

to be met, it makes sense because the elements in both sides are central. 
With all this in place, the following relation 

— ^n— 1 — Q UnQ 1 

with 



Q = hb 2 . ■ ■ 6iv_i b n = 



proves to hold for all n E Z. So, with automorphism" : u n \— > n n _i inter- 
preted as a jump in time, the above relation replaces the (lower) Volterra 
system with a free and fully discrete Heisenberg quantum motion, employing 
the element Q as an "evolution operator". 

1.6. Proof. We just honestly start from n = 1: 

00 

Q- 1 u 1 Q = Q- 1 u 1 b 1 b 2 b 3 ...b N _ 1 =Q-\ q a2+2a u a 2 u 1 b 3 ...b N _ 1 

a=— 00 

= Q~ l b l b 2 q ^2 1 U\b% ■ ■ ■ b N _x = Q~ 1 b 1 b 2 b 3 u 2 ~ 1 u 1 u 3 b 4: . . . b N _ 1 = ■ ■ ■ 

- every step adds another factor, in the end parity of N and that center- 
reducing condition (||) make all the difference — 

• • • = Q^Qi^u^ . . . u Af _ 2 ) _1 u 1 u 3 . . . u N _ 1 = u N . 

We may seem to have done only one very particular case, but, luckily, the 
rest is little more than tautology. It immediately follows that Q~ l b\Q = bjy, 
which can be written as 

Q = b^Qbx = b 2 b 3 ...b N , 

which in turn just means that Q = Q and therefore anything good for n = 1 
at once becomes just as good for all n. That is it. 

In the process, a seeming contradiction between the "open-ended" ap- 
pearance of the evolution operator and the cyclic nature of its action has 
resolved itself. Since 

Q = hb 2 . . . b N _ 1 = Q = b 2 b 3 . . . b N = Q = ■ ■ ■ = b N b N+ i . . . b N - 2 , 
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whichever of these forms one chooses to express Q in, it always remains an 
ordered product but each time starts from another point, in this sense it 
does not depend on the starting point, or rather does not have one. The 
next subsection offers an explanation of this miracle. 

1.7. It has to do with braids. Let us compile a list of what we know 
about the Vs. They are periodic, b n+ N = b n , they are this "local" 

b m b n = b n b m if \m - n\ ^ 1 (mod N), 

and they satisfy "global" relations which closed the previous subsection 

&l6 2 • • • &JV-1 = &2&3 • • • 6jv = • • • = &tv6tv + i . . . &AT-2 ■ 

Let us now consider these as defining relations and identify the emerging 
group. It is Bn, the group of braids of TV strings in 3 dimensions, alterna- 
tively defined by an exhaustive list of (Artin's) relations 

b n b n -ib n = 6 n „i6 n 6 n _i b n b m = b m b n if |n - m\ ^ 1 

imposed on N — 1 generators 6i, 62, - - - , b^-i- We leave it as an (instructive) 
exercise to check that our list is equivalent to Artin's, provided the leftmost 
of our "global" relations is read as 

b N = (b 2 . . . b N - 1 )~ 1 b 1 b 2 . . . &jv-i 

and understood as the definition of 6 at. Just in case, let us warn against 
mistaking bjy for the iV'th generator of B^ + \, ours is just an element of 
-Bjv explicitly defined above. Those familiar with the braid group must have 
already recognized in it the first string crossing the very last one behind all 
the others, let us picture it for N = 4: 



XII -IX I -I IX -X) 



A picture of the evolution operator 

Q = &1&2&3 = &2&3^4 = hb4h = &4&1&2 



now explains better than words how it manages to be ordered and cyclic at 
the same time. 

1.8. Lattice Virasoro algebra. Now we look for a generalization of this 
construction by natural interpretation of the "raising". So, we want the 
lattice Virasoro algebra actually to be a group. If that is to be, the "global" 
relations ought to remain the same as in the braid group 

did,2 ■ ■ ■ dN-i = d 2 ds . . . = ■ ■ ■ = d]ydN+i ■ ■ ■ d 2 N-2 , 
or equivalently 

(|") d n = d n -i = Q~ 1 d n Q Q = dxd.2 ■ . . d N -x , 

but this time each generator should interfere not only with the nearest neigh- 
bours but also with the second-nearest ones 

d m d n = d n d m if |m — n\ ^ 1 and 2 (mod N), 
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just like it was in bracket (j5[) of course. This by itself implies 

d n dn-2dn-id n dn+i = d n -2dn-id n dn+id n -i 

instead of Artin's relations, but we feel that these are not tight enough, so 
we voluntarily split each of them in two 

d n dn-2d n -id n d n+ i = d n -2d n dn-ldn+l = d n -2d n -i 

dn dn + 1 d n — i , 

to end up with weird relations from 47] 

d n +id n -id n d n +i = d n -id n +id n 

d n +i — d n -\d n d n+ id n ^i. 



(I") 



This coup de force will find its justification in the treatment of the Liouville 
model below. 



2. Difference-difference Liouville equation 

2.1. Liouville formula. For aesthetical reason alone, its difference approx- 
imation can only be 



2 20 {g> - g)(f3> - (3) 
~ {a' - (3'){a - f3) ' 



with primes now denoting finite shifts of arguments by A. Appearance aside, 
we have already mentioned that the Liouville formula manifests invariance of 
<f> under simultaneous point-wise Mobius transform of chiral halves a and (3. 
The so called cross-ratio employed in our difference scheme has just the same 
property. 

To make of this a lattice formula we draw a square j, k lattice and put 




Xjk 



ft 



|(i+fc+i) 



(j+k-l) 



a hU-k+i) 



%j+fc-l)) 



in the vertices with j + k = 1 (mod 2), like this: 

k p 3 ai fa a 2 




So, the x' s occupy sites marked by bullets, the empty circlets will this time 
remain empty. The second /c-axis at j = 8 reminds that the lattice covers a 
cylinder rather than a plane, Xj+2N,k = Xjk , with not necessarily equal 
four of course. 
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2.2. Discrete Liouville equation. It not only ideologically justifies the 
chosen discretization but also helps make the following action more en- 
tertaining. In continuum, the Liouville equation may be considered as a 
compatibility condition for the Liouville formula overloaded with Mobius- 
invariance. It ought to be just the same on the lattice, so let us find out 
what "equation of motion" the x's might solve. We pick four of them next to 
each other to figure out that they are four cross-ratios made from the total 
of six a's and /3's. That is one too many: 4 + 3 symmetry parameters 
- 3 a's - 3 /3's = 1. The missing one turns out to be 

(7) Xi,k + lXj,k-i = jz - ^"'ffi + 'v T j + k = (mod 2) 

I 1 + Xj-l,k){± + Xj+l,k) 

which so becomes our favourite lattice Liouville equation — replacing Hi- 
rota's original 

n-j,k+i n j,k-i — , , r 7 • 

The two are connected by a simple change of variables 

Xjk = hj-i^hj+ifi, 

which quarters the /i's in the empty sites of our lattice of course. However, 
that change would also split cross-ratios, which we want to avoid, at least 
in this paper. 

The stress-energy connection (FJ) also lends itself to a cross-ratio treat- 



ment. The Schwarz derivative inevitably becomes another cross-ratio [ 1 7 [ 

A K+2-Qti+l)K-"n-l) 

(«n+2 - a n -i)[a n+ i - a n -i) 
then another counting argument leads to 

(8) — ^— = ( 1 + Xj ,k-l + %*=±) ( 1 + X W + M±1A - 



u \(j-k) V Xj-l,k/ V Xj+2,fc-i, 

2.3. Poisson bracket. Two rows of x's make perfect Cauchy data 

X2n = X2n,-1 X2n+1 = X2n+1,Q ; 

just in case one may directly check that the lattice canonical bracket of [21] 
{X2n±l,X2n} = X2n±lX2n { X j,Xi} = if \j - i\ # 1 (mod N) 

both a) reproduces itself as the x' s evolve according to the lattice Liouville 
equation (0) and b) translates, by means of formula (|8|) at k = 0, into the 
lattice Magri bracket (|5f). 

3. Formal quantization 

3.1. Algebra of observables. In a virtual replay of subsection 1.5, we in- 
troduce a formal quantization constant q, replace the above lattice canonical 
bracket by Weyl-style exchange relations 

X2n±lX2n = <? 2 X2nX2n±l XjXi = XiXj ^ \j ~ *| ^ 1 (mod 2N), 

opt for even N and impose an additional condition on central elements 
(9) 

X1X3 X5X7 • • • X2JV-3X2AT-1 — X2X4 X6X8 • • • X2iV-2X2Af • 
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3.2. Quantum lattice Liouville equation. We can only afford a Heisen- 
berg evolution of observables, so let us again consider Xj as initial data 
(X2n,-l = X2n, X2n+i,o = X2n+i) and determine elements Xj,k+1, j + k = 
(mod 2) step by step using a slightly "quantized" lattice Liouville equation 

/m\ Q 2 Xj-l,kXj+l,k 

( W > Xj,k+lXj,k-l = -rr— — 77--. r • 

l 1 + qxj-i,k){t + qxj+i,k) 

We deliberately present the r.h.s. as a ratio to stress that all the factors 
there commute with each other, on the contrary, those in the l.h.s commute 
neither with each other nor with the r.h.s.. 

3.3. Shift and evolution operators. If we guessed the equation right, 
there must exist the "evolution operator" , that is an element K such that 

Xj,k+l = K^Xj^-iK. 

Needless to say, these relations will hold everywhere if and only if they do 
so for k equal and —1 where they become 

js q 2 X2n-2X2n 
X2n~lKx2n-lK ~ 



(1 + <?X2n-2)(l + QX2n) 

is q 2 X2n-\X2n+l 
K X2nKx2r, ~ 



(l + gx2n-i)(l + qx2n+i) 



Following [19, 21] we shall produce that element, with explicitly singled out 



d'Alembert part good for 

X2n-lK OQ X2n-lK^ = q 2 X2n-2X2n K^X2nK 00 X2n = q 2 X2n-lX2n+l, 



and complete with the "shift operator" J moving the x' s m the spatial 
direction 

Xj+l,k = J 1 Xj-i,kJ- 
Here follow explicit formulas for these shift and evolution operators 



K OD = VU K = E 2 KE 1 J = VU- 1 



where 



N N 

E l = ll <X2n-l) E 2 = J] e( X 2n) 

n=l n=l 

U = 0(qXi 1 X2)0(qX3 1 X4) ■ ■ ■ 6(qX 2 N-3X2N-2) 

y = G(qx2N-ix 2 N-2) e (qx2N-3X 2 N^) ■ ■ ■ d (qx 3 x 2 1 )- 

Of course, U and V have everything to do with braids but we will not go 
into that. The two special functions involved 



e(z) = f[(l + q 2 P +1 z) = (-qz; q 2 ) c 



p=0 



{z) =e{z)e{z~ l ) = const- q l? z p 
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are quite special indeed but all we want to know about them, for now at 
least, are the beautiful functional equations 

e(gz) 1 6{qz) 1 

e(q- 1 z) ~ 1 + z Oiq^z) ~ ~z 

fulfiled by the former, and the ensuing equation on the latter which has 
already been used, albeit implicitly, in subsection 2.4. This time those func- 
tional equations and condition (|9|) gradually translate into relations 

E l l X2n-l E l = X2n-1 E \ l X2n E l = X2n( 1 + QX2n-l)( 1 + <7X 2 n+l) 

E 2 X X2n E 2 = X2n ^XSn-l^ = ((1 + WC&l^X 1 + <lX2n))~ X X2n-1 

U~ 1 X2nU = Xan-1 U~ l X2n+lU = ^ X2n-lX2n+lX2n 

V ~ 1 X2n V = X2n+1 V ^X2n-l V = f X2n-lX2n+lX2n 

which combined make J and K satisfy what they have to. We omit the 
calculation but a remark is in order. 



3.4. Odd N or no condition (|9|). In these cases it all fails, and subsec- 
tion 1.6 gives an idea why. In fact, they are better served by those cross- 
ratio-splitting variables h and Hirota's lattice Liouville equation which we 
were so quick to discard. We are planning to return to this issue elsewhere. 



3.5. Chiral evolution operator. We define it as 

Q = UEx , 

and it indeed moves the x' s i n the right direction 

Xj,k+i = Q^ 1 Xj+i,kQ 

and equals V J _1 K, in the sense that Q 2 = J~ X K. So, the total momen- 
tum P, Hamiltonian H and the chiral Hamiltonian \{H — P) of the original 
equation have finally become shift-evolution-operators J, K and Q of the 
quantum lattice equation, but we have yet to find out in what sense, if at all, 
this Q coincides with that hypothetical quantum lattice group-like counter- 



part of Lq which we called Q in section L8 . The answer does not come easy 
but in the end few carefully placed q's do it again, and so does the magic 
function e introduced in subsection |3.3|. It turns out that 



(11) Q = d\d,2 ■ ■ ■ d N -i, 

where Q is UE\ and the d's are 

dn = e(^ _1 (l + qX2n + X2nX2n-l)( 1 + VX2n+l + X 2 n+2X2n+l) ~ ^) ■ 

As the notation suggests, these also satisfy relations (Q 1 ') and (|f). It is 
instructive now to see, that argument of function e in this formula is a 
natural quantisation of expression (|8|). 
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3.6. Proof of eqn (pTl|) . First, we recall the two identities satisfied by 
e- function: 

(12) e(u)e(v) = e(u + v), e(v )e(u) = e(v + u + qvu), uv = q 2 vu. 
Next, we have the identity 

e(Xi(l + qX2 1 ))<3 = Q<X2N^ + ^-i)" 1 ) 
which is a consequence of the relations in subsection |3.3] satisfied by U and 
E± operators. It is also easily checked that for any j 

0(qX2j~-iX2j)e(x2j-i) = e(xjy_i(l + %/))e((l + qx%-i)X2j)- 
Now, using these relations, we have 

Q = (e(xi(l + q X2 1 ))r 1 UE 1 e(x2N^ + QXm-l)' 1 ) 

= (e( X i(l + qX2 1 ))r 1 I II 8(<lXy-iX2jHX2j-i) 

\l<j<N 

x e(x2N~i)e(x2N( 1 + QX2N-1)' 1 ) 
= (e(xi(l + qX 2 1 ))r 1 [ II zfaj-iil + qX^Mil + qXy-dXij 

\l<j<N 

X e(9X2Ar-lX 2 Jv) e (X2Ar-l) 

= II e (( 1 + ^2/-l)X2j)e(X2j+l(l+gX 2 /+2)) = II 
l<j<N l<j<N 



thus obtaining eqn (|TT|) . 

The proof of relations (g") is done in [47] and is also based on identities 
©■ 

4. DUALIZATION 

4.1. Involution. The formal quantization of section ^ can be put in the 
usual framework of quantum mechanics if we assume that both the formal 
parameter q and the generators x are exponentials 

q = e mT Xj = e - 2 ^i, 

and the new generators (fj have the commutation relations 

[f2n±l, P2n] = [fjiVi] = if |j - i| ^ 1 (mod N) 

independent of r, and so could be taken as selfadjoint operators 

V>] = V>y 

The rest depends on r of course. If the established formula 

c=l + 6(r + i + 2) 

(relating the coupling constant of the continuous Liouville field theory to the 
central charge of the corresponding representation of the Virasoro algebra) 
has something to do with our lattice theory, then three cases have to be 
considered: a) c < 1 «-> r < 0; b) c > 25 <-> r > 0; and c) 1 < c < 25 <-> 
It I = 1. The first two, although lead to pretty normal reality conditions 
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Xj = Xj 1 an d Xj = Xj respectively, also put q on the unit circle, which is 
the last thing we want right now. That leaves c). 

4.2. Change of function e. So, |r| = 1, also, since c does not distinguish 



between r and 1/r, let for definiteness 3r > 0. Following [13|, consider the 
function 

(13) /(*) = e^iz/y/r) = ) * ' J°° , 

which uses our former favourite e as the numerator but divides it by itself 
but with suitably altered arguments. In Appendix we collect some of the 
properties of this function. Here we remark in particular that f(z) satisfies 
the same functional equation 

(u) *Lti) = 1 

1 j /(z-§) l + e 2 «" 

as e(e 2mz ) did, but this time we also have 

|/(*)l = l 

on the line z = rz, simply because the numerator and denominator of / are 
complex conjugate of each other on that line. There is an easy profit to be 
had from that. 

4.3. Unitarity. Replace e by / in every factor of every shift-evolution op- 
erator, all those little factors and big operators will at once become unitary 
but all the relations we have found them to satisfy will remain intact, for 
they rely only on that one and only functional equation. For instance, let 
us write out the so upgraded chiral evolution operator: 

N 

Q = K 2(jV-l) e «(vn-V>2) 2 e 7Ti( V 3-^4) 2 _ _ _ ^i(V2N-3-V2N-2) 2 JJ /(fr/F^-l), 

71=1 

where it is already taken into account that 

f(z)f(-z) = K 2 e-^\ 

with k = /(0) of course. 

Mission accomplished, our lattice Liouville model has finally turned from 
an algebraic fantasy into a quite material unitary theory. Moreover, we shall 
momentarily see that we get a new feature. 

4.4. Dual Liouville equation. Let us permute the factors in the l.h.s of 
the lattice Liouville equation ( |l~0| ) 

Xj,k-lXj,k+l = Xj,k-1 [Xj,k+lXj,k-l) Xj,k-l 



Q 2 Xj-l,kXj+l,k \ -l 

Xhk ~ X V(l + QXj-i,k)(l + qXj+l,k)J Xhk - 1 

Q~ 2 Xj-l,kXj+l,k 

v 1 + 9 _1 Xi-i,fe)(l + Q^Xj+hk) 
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and treat the both sides with Hermitean conjugation, remembering that 
= Q^ 1 of course. Since 

Q = q 1 = e -r = X j = Xj = e ^ = x] , 

the resulting equation reads 

~2 ~ 

m) = (1+ 4^aT4 + i, fc )- 

It is plain to see that it is the same equation but with q and x' s replaced 
by q and %'s, this is called duality. Indeed, instead of conjugating things, 
we might instead explore the fact that those x's satisfy the same relations 

X2n ± lX2n = ?X*nX»±l 

as the original do, except with q instead of q. Moreover, the two sets 
commute with each other 

XjXi = XiXj ■ 

The algebras they generate (call them A T and Ai) form two factors in the 

r 

algebra B generated by the tp's, and leave no free space, in the sense that 
B = A T Ai , sort of. Either way, such a bisection is well served by the 

T 

function /, which fittingly satisfies the dual functional equation 

(0 ) T7 TT = 2^77 

1 + e— z 

to complement the original one (0). Equation (0*) can now be derived 
in exactly the same way as (|i~0| ) was, and regarded as not just a conjugate 
clown of the latter but as an equal dual equation. 

4.5. Baxter equation. It reads 

t T (X)Q(X) = Q(X + \) + (e^ iTX + l) N Q(X - ±), 

where t T (X) and Q(X) are two families of elements (of A T and B respectively), 
which all commute with each other 

[t T (X),U(»)} = [Q(X),tr{fi)] = [Q(A),Q( M )]. 

In this paper we do not define those families explicitly and do not derive 
the equation. Let us only mention that the chiral evolution operator Q is in 
fact Q(X) evaluated at a particular value of A. The duality symmetry of our 
construction implies that there is another dual equation having the form 

U (A)Q(A) = Q(X + + (e imX + if Q{X - §). 

T 

Thus, the Baxter equation gets accociated with the modular lattice. 

We stop here and postpone the study of Baxter equations to the next 
paper of this series. 
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5. Conclusions and comments 

We have shown, that the appropriate lattice regularization of the quan- 
tum Liouville model allows to enter the "forbidden region" of the Virasoro 
central charge 1 < c < 25. The key for this is using the double family of 
dual quantum fields, which are not selfadjoint, but normal and adjoint to 
each other. The results are still rather modest, but proving the point quite 
persuasively. 

The real problem is that of findng the spectrum of the model. One viable 
approach to this can be based on the Baxter equation. In the next paper 
of this series we plan to deal with this equation. The idea of its derivation 
originally goes to Baxter himself . In j|, |4| it was realized in the context of 
the chiral Potts model, and in ]|, [?J in the context of quantum conformal 
field theory. The dual Baxter equations in a different context were also 
considered recently by Smirnov f39fl . 

In parallel construction for the WZNW model |]12[ it was shown how to 
separate the zero mode problem which reduces the problem of finding the 
spectrum of the primary states to a problem with finite number of degrees 
of freedom. It is possible that similar construction can be found here. If 
so, one will be able to bypass Baxter equations for the investigation of the 
spectrum. 

Liouville model is known to be a contraction of the massive Sine-(Sinh-) 
Cordon model to which most of our considerations are also applicable J2(]]. 

In a series of papers H || 0] the quantum KdV equation was considered 
without recourse to lattice. As the Liouville and KdV models are close 
relatives, it is no wonder that one can see many similarities between those 
papers and our text. Moreover, the discrete variant of quantum KdV was 



already considered in |46|| . We believe, that using lattice allows to use the 



duality in full strength and, in particular, go to the strong coupling region. 



In papers [27, 28 1 arguments were raised for the exceptional values of the 
central charge c = 7, 13, 19. From our point of view, these values are distin- 
guished only by giving our modular lattice supplementary symmetries (r cor- 
responds to elliptic fixed points of the modular figure: e l7r / 3 , e l7f//2 , e j27r / 3 ). 
There is no reason for us to exclude other values of r on the unit half-circle. 

Finally, we cannot help stating that for us the function f(z) defined in 
eqn (|13|) seems to be the real cornerstone of the theory of quantum integrable 
models. It appears in basic objects of the dynamical theory: the evolution 
and the Baxter operators. Its close relative as a function of rapidity defines 
the Zamolodchikov's factorized S-matrices. It is also indispensable in the 
theory of form- factors We believe that the full content of duality, 

hidden in this function, is still not explored to full extent. In the Appendix 
below we describe some of its remarkable properties. 



6. Appendix: the non-compact quantum dilogarithm 

Let complex b have a nonzero real part !R& ^ 0. The non-compact QDL, 

et„(z), z G C, < 



c^iib + b- 1 )^, 



STRONGLY COUPLED QUANTUM DISCRETE LIOUVILLE THEORY 15 

is defined by the formula 

(\ f +oc e~ 2izx dx \ 

( 15 ) eb ^ ~ ex P 7 / ■ u ( u\ ■ \>t iu\ ' 

\4 J_ OQ smn(xo) smh(x/6)x / 

where the singularity at x = is put below the contour of integration. 
This definition implies that e b (z) is unchanged under substitutions b — > 
b — > —b. Using this symmetry, we choose b to lay in the first quadrant of 
the complex plane, namely 

m > 0, 96 > 0, 

which implies that 

9r > 0, t = b 2 . 
This function has the following properties. 

6.1. Functional relations. Function ( |i~5| ) satisfies the 'inversion' relation 

(16) e b {z)e b (-z) = e ^ 2 -^+H)/^ 
and a pair of functional equations 

(17) e b (z - ib ±l /2) = (1 + e 2 * zb±1 )e b (z + ib^/2). 

The latter equations enable us to extend the definition of the QDL to the 
entire complex plane. 

When b is real or a pure phase, function e b (z) is unitary in the sense that 



(18) e b (z) = l/e b (z). 

If selfadjoint operators P and X in L 2 (W) satisfy the Heisenberg commutation 
relations 

(19) [P,X] = ^, 
the following operator five term identity holds: 

(20) e b (P)e b (X) = e b (X)e b (P + X)e 6 (P). 

For real b this can be proved in the C*-algebraic framework Q. We will prove 
it for complex b. The case of real b then will follow by continuity. 

6.2. Analytic properties. We can perform the integration in (15) by the 
residue method. The result can be written as ratio of two q-exponentials 

(21) e b (z) = (e 2 -( 2+c ") fe ;g 2 ) 00 /(e 2 -( 2 - C6 ) fe " 1 ;Q 2 ) 00 , 
where 



„27r6 2 - „— i-wb 2 



q=e , q = e 

thus reprpoducing definition (^) in the text. Formula (21) defines a mero- 
morphic function on the entire complex plane, satisfying functional equa- 
tions (|l^) and (0), with essential singularity at infinity. So, it is the ana- 
lytical continuation of definition (|l5|) to the entire complex plane. It is easy 
to read off location of its poles and zeroes: 

zeroes of {e^z))^ 1 = {t(q> + mib + nib -1 ) : m,n S Z>q}. 



S.L. Woronowicz: private communication, 1998 
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The behavior at infinity depends on the direction along which the limit is 
taken: 



(22) e b (z) 



where 



\z\— >oo 



1 I arg(»| > f + arg(fe); 

e inz 2 -in(l+2ci)/6 | aTg ( z )\ < | _ arg (fr) ; 

(<f ;<f WB^" 1 -?; -6~ 2 ) |argz-7r/2| <arg&; 

k ©(ibzib 2 )/^ 2 -^ 2 )^ |argz + 7r/2| < arg6, 



0(z; t) = Y^ e inTn2+2ninz , > 0. 



Thus, for complex b, double quasi-periodic ^-functions, generators of the 
field of meromorphic functions on complex tori, describe the asymptotic 
behavior of the non-compact QDL. 

6.3. Integral Ramanujan identity. Consider the following Fourier inte- 
gral: 

(23) ^(u,v,w)= [ eb ^ + U ] e 2mwx dx, 



e b {x + v) 
where 

(24) Q(v + c b )>0, %(-u + c b )>0, 3f(u - u) < <3u> < 0. 



Restrictions (|24j) actually can be considerably relaxed by deforming the in- 
tegration path in the complex x plane, keeping the asymptotic directions of 
the two ends within the sectors ±(| argx| — 7r/2) > argfr. So, the enlarged 
in this way domain for the variables u, v, w has the form: 

(25) | arg(iz)| < ir — arg6, z £ {w, v — u — u>, u — v — 2c b }. 

Regarding e b (z) as a 'non-compact' analogue of the q-exponent (x;g)oo, 
definition ( p3[) can be interpreted as the corresponding integral counterpart 
of the Ramanujan sum: 



Hpi(x,y,z) = 



[x] q)n ^ n 



The latter is known to be evaluated explicitly, the result being the famous 
Ramanujan summation formula: 

, / \ (q;q)oo(y/x;q) 00 (xz;q) 00 (q/xz;q) 00 
m{x,y,z) - 



(y; q)oo(q/x; q)oc(z] q) 00 (y/xz; q) c 



Remarkably, integral fl23| ) can be evaluated explicitly as well. Indeed, using 
the residue method, we easily come to the following result: 

(26) <S(U, V, W) = ~ V ~ Cb)eb{ - W + Cb) c -2^^+ C ,)+^(l-4cg)/12 

e b (u - v + w - c b ) 

( 27 ) = e b(v -U-W + C b ) e -2 7 rtt t »(M-c i ,)-t7r(l-4^)/12 

e b (v -u + c b )e b (-w - c b ) 
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where the two expressions in the right hand side are related to each other 
through the inversion relation (|i~6|). The similarity of this result with Ra- 
manujan sum becomes very transparent if we rewrite the latter in the form: 

(os\ y (£/g n ;g)oo » _ (y/x;g)oo(q/z;q) 00 e q (xz) 2 
1 ' ^N";g)oo (y/xz-^u 9 q (x)e q (z) [q ' q) ^ 

where the ^-function is defined by 

(29) 6 q {x) = (g;g)oo(x;g)oo(^;g)oo = ^g" ( "- 1)/2 (-xr. 

Comparing the inversion relation (16) with eqn (|29|), we conclude that the 
non-compact analogue of the ^-function is the Gaussian exponent, and the 
structures of eqns (pg|) and (pB"|) are now quite similar. 



6.4. Fourier transformation of the QDL. Certain specializations of 
^f(u, v, w) lead to the following Fourier transformation formulas for the QDL: 



(30) Mw)= / e b {x)e 2mwx dx = y{0,v,w)\ v ^ oo 
Jr 



and 

-1 2-iriwx 



(31) = / (e b {x)Y l e Zmwx dx = ^{u, 0, w)]^^ 

Jr 

= e -2~+-(l-4c?)/12 e?)(u; + Cb) = e in W *-iTr(l-4cl)/l2 /eb (_ w _ ^ 
The corresponding inverse transformations read: 

(32) (e fe (x)) ±1 = / Mv)e- 2mxy dy, 

Jr 

where the pole at y = is surrounded from below. 

6.5. Proof of the Pentagon identity. Using formula (32) and commuta- 



tion relation (19), we equate the coefficients of the operator terms 



in the pentagon relation (|20|), the result being an integral identity: 



where the singularities at z = x, z = y are put below, and at z = 0, 
above the integration path. Now, multiplying both sides of this identity by 
exp(— 2myu), integrating over y, and using (^), we obtain 

(j) + {x)e b {u-x) = e b {u) I <f) + {z)4> + {x - z)e i7Tz2 - 2iriuz dz. 



Using (pOQ, we rewrite it equivalently 

(33) e b (u-x) ^(i-4cg)/i2 
e b (-x - c b )e b (u) 



e b (z + c fe ) 
e b (z - x - c b ) 



e -2viz(u+c b ) dz = _ x _ c&j _ u _ 
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which is a particular case of 
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